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ABSTRACT 

Nucleus-encoded ribonucleases and RNA-binding 
proteins influence chloroplast gene expression 
through their roles in RNA maturation and stability. 
One mechanism for mRNA 5' end maturation posits 
that sequence-specific pentatricopeptide repeat 
(PPR) proteins define termini by blocking the 5 '-*3' 
exonucleolytic activity of ribonuclease J (RNase J). 
To test this hypothesis in vivo, virus-induced gene 
silencing was used to reduce the expression of 
three PPR proteins and RNase J, both individually 
and jointly, in Nicotiana benthamiana. In accord- 
ance with the stability-conferring function of the 
PPR proteins PPR10, HCF152 and MRU, accumula- 
tion of the cognate RNA species atpH, petB and 
rbcL was reduced when the PPR-encoding genes 
were silenced. In contrast, RNase J reduction 
alone or combined with PPR deficiency resulted in 
reduced abundance of polycistronic precursor tran- 
scripts and mature counterparts, which were 
replaced by intermediately sized species with het- 
erogeneous 5' ends. We conclude that RNase J de- 
ficiency can partially mask the absence of PPR 
proteins, and that RNase J is capable of processing 
chloroplast mRNAs up to PPR protein-binding sites. 
These findings support the hypothesis that RNase J 
is the major ribonuclease responsible for maturing 
chloroplast mRNA 5 r termini, with RNA-binding 
proteins acting as barriers to its activity. 

INTRODUCTION 

Chloroplasts arose from the phagocytotic uptake of a 
photosynthetic bacterium by a primitive heterotrophic 
eukaryote. Following endosymbiosis, the bacterium sur- 
rendered its autonomy through massive gene transfer to 
the host nucleus (1). Although the chloroplast retains its 



own transcription and translation apparatus, gene regula- 
tion is predominately conferred by nucleus-encoded 
proteins targeted to the chloroplast (2). The majority of 
chloroplast-encoded genes are arranged into operons or 
clusters, and the resulting polycistronic transcripts are fre- 
quently processed into mono- and oligocistronic RNAs. 
These processing events involve a mosaic of trans-factors, 
namely endo- and exoribonucleases of mainly prokaryotic 
ancestry together with RNA-binding proteins of more 
recent origins (3). Here, the cooperative functions of ribo- 
nuclease J (RNase J) and the pentatricopeptide repeat 
(PPR) family of RNA-binding proteins were investigated. 

RNase J was discovered in Bacillus subtilis, where it was 
demonstrated to be an essential protein harboring both 
endonuclease and 5'-*3' exonuclease activities (4,5). It is 
found in Archaea, Gram-positive bacteria, cyanobacteria 
and chloroplasts but is absent from Escherichia coli. 
Nuclear RNase J-encoding genes are ubiquitous in 
plants, and Arabidopsis RNase J is targeted to the 
chloroplast, possesses dual ribonuclease activities (6) and 
is essential for embryo development (7). Tobacco 
{Nicotiana benthamiana) plants made deficient in RNase 
J by virus-induced gene silencing (VIGS) exhibit a remark- 
able accumulation of antisense RNAs (asRNAs), which 
duplex with sense transcripts and impair their translation 
(6). These asRNA substrates of RNase J are hypothesized 
to originate from inefficient chloroplast transcriptional 
termination. RNase J, therefore, participates in a chloro- 
plast RNA quality-control mechanism by eliminating 
read- through transcription products. 

Recent studies have also implicated RNase J in a 
chloroplast mRNA maturation pathway directed by 
PPR proteins. PPR proteins are characterized by helical 
repeat structures formed from tandem arrays of a degen- 
erate 35 amino acid motif and are prevalent among land 
plants but absent from bacteria (8). Many are organelle 
targeted, bind RNA in a gene-specific manner through an 
amino acid code (9,10) and confer transcript stability by 
obstructing ribonucleolytic digestion and defining mature 
mRNA termini (11). In the case of polycistronic mRNAs, 
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PPR proteins bind overlapping transcripts relieved from 
precursors by multiple endonucleolytic cleavages within 
intergenic regions and block exonucleolytic processing 
from both the 3' and 5' directions (12,13). PPR 10 and 
HCF152 bind the atpI-atpH and psbH-petB intergenic 
regions, respectively, and shield transcript 5' and 3' 
termini mapping to those regions (13,14). PPR proteins 
are also known to define monocistronic transcript ends 

(15) . MRL1, for instance, binds a sequence in the 5' 
untranslated region (UTR) of the primary rbcL transcript, 
leading to the accumulation of a shorter processed form 

(16) . The enzyme responsible for maturing PPR protein- 
defined 5' ends has yet to be determined. As a commercial 

exonuclease can carry out this function in vitro (17), 
RNase J is a prime candidate for this role in vivo as the 
only identified chloroplast 5 / ^3 / exoribonuclease, as pre- 
viously speculated (13). 

To evaluate the involvement of RNase J in PPR 
protein-mediated 5' mRNA maturation, genes encoding 
PPR10, HCF152, MRL1 and RNase J were silenced, 
independently and in pairwise fashion, by VIGS in 
N. benthamiana. Deficiency of single PPR proteins 
resulted in diminished accumulation of processed 
cognate RNAs, as expected, an effect that could be par- 
tially masked by deficiency of RNase J. Based on these 
results, we conclude that chloroplast RNase J trims 
mRNA 5' termini to mature forms defined by bound 
PPR proteins. 

MATERIALS AND METHODS 

Plant material 

Nicotiana benthamiana plants were grown in a greenhouse 
at 25°C under a 14-h light/ 10-h dark regime. 

Design of VIGS target regions 

Arabidopsis genes for PPR10 (At2gl8940), HCF152 
(At3g09650) and MRL1 (At4g34830) were submitted 
as BLAST queries of the Nicotiana tabacum Unigene 
sequence set of the Sol Genomics Network (18). 
The largest open reading frame of each Unigene result 
was translated and aligned with the corresponding 
Arabidopsis PPR protein sequence using ClustalW2 (19), 
verifying that the orthologous N. tabacum genes had been 
identified. For VIGS target regions, two unique stretches 
of 250-350 nt for each N. tabacum sequence were selected, 
which returned only a single hit on submission to siRNA 
Scan (http://bioinfo2.noble.org/RNAiScan.htm). All 
primer sequences are given in Supplementary Table SI. 

Plasmid construction and VIGS infection 

One or two target regions for PPR10, HCF152, MRU 
and RNJ were PCR-amplified from N. benthamiana 
DNA. The amplicons were cloned into pCR®8/GW 
(Invitrogen) as both single inserts and RNase J-PPR 
protein tandem inserts. The VIGS targets were then 
inserted into pYL170/GW (20) and transformed into 
Agrobacterium tumefaciens GV3101 by electroporation. 
As a control, the empty vector pYL170, along with the 



expression component of the bipartite VIGS system, 
pTRVl, were introduced into Agrobacterium. VIGS was 
carried out as previously described (6). Given that genes 
silenced with either of two targeting regions consistently 
yielded indistinguishable molecular phenotypes results 
from only one VIGS targeting region per gene or gene 
combination are presented here. 

RNA extraction, 5' rapid amplification of cDNA ends and 
quantitative RT-PCR 

Leaf tissue was harvested 18 days after VIGS infections, 
and RNA was isolated with TRI reagent (Sigma-Aldrich) 
per the manufacturer's instructions. Three micrograms of 
RNA were treated with 2 U DNase I (Ambion) for 30 min 
at 37°C. The reaction was terminated by phenohchloro- 
form extraction, followed by ethanol precipitation. 
Treatment with tobacco acid pyrophosphatase (TAP) 
(Epicentre; 0.5 U per 3 ug RNA) was carried out at 37°C 
for 1 h with 40 U RNaseOUT (Invitrogen). The RNA was 
re-purified and incubated with 1.7 ul of RNA adapter 
(10 uM; see Supplementary Table SI) at 65°C for 5 min, 
briefly chilled on ice and then added to the ligation 
mixture [ImM ATP, 4U RNaseOUT, lx ligase buffer, 
0.5 U T4 RNA ligase (New England Biolabs)]. 
Superscript III (Invitrogen) was used for first-strand syn- 
thesis, primed with random hexamers, in accordance with 
the manufacturer's instructions. Reverse transcription 
products were amplified by PCR with 1 U of GoTaq 
(Promega) using primers specific for the RNA adapter 
and the gene of interest. The following conditions were 
used: 2 min at 95°C; 30s at 95°C, 30s at 54°C, 70s at 
72°C, repeated for 30 cycles; and 2 min at 72°C. Rapid 
amplification of cDNA ends (RACE) products were 
separated in 1.5-3% agarose gels, and DNA fragments 
were gel purified and ligated into pGEM-T Easy 
(Promega) for sequencing. For quantitative RT-PCR 
(qRT-PCR), strand-specific cDNA synthesis was per- 
formed using the 3' qPCR gene-specific primers 
(Supplementary Table SI), and PCR was conducted as pre- 
viously described (21) using the same reference genes, 18S 
rRNA and glyceraldehyde 3-phosphate dehydrogenase 
(GAPDH). Standard errors were calculated from three bio- 
logical replicates. 

RNA gel blot analysis 

Total RNA (0.5 ug/lane) was fractionated in 1.2% formal- 
dehyde-agarose gels and transferred to Amersham 
Hybond-N+ membranes (GE Healthcare). Strand- and 
gene-specific probe templates were PCR-amplified from 
N. benthamiana DNA (see Supplementary Table SI) and 
gel purified. RNA probes were generated by T7 RNA 
polymerase (Promega), in accordance with the manufac- 
turer's protocol, with 40 uCi of oc 32P -uridine-5 / - 
triphosphate (UTP). RNA gel blots were prehybridized 
as previously described (6), hybridized with riboprobes 
for 14h at 68°C, washed twice in lx SSC, 0.6% SDS 
for lOmin at 68°C and twice with 0.1 x SSC, 0.6% SDS 
for 45 min at 68°C. DNA-probed RNA gel blots were per- 
formed as previously described (22). 
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RESULTS 

RNase J processes atpH 5 ends defined by PPR10 

To study the function of RNase J in vivo, expression of the 
enzyme was repressed in N. benthamiana by VIGS using a 
bipartite tobacco rattle virus system (23). As we previously 
reported (6), infected plants deficient for RNase J 
exhibited a chlorotic phenotype (Figure 1A), indicative 
of impaired photosynthetic capabilities. Indistinguishable 
physiological and molecular phenotypes were seen when 
VIGS targeted either of two unique regions within the 
RNase J-encoding transcript, indicating minimal off- 
target silencing. Efforts to use immunoblot analysis to 
quantify RNase J deficiency were unsuccessful; however, 
our earlier work showed that VIGS reduces the transcript 
level for RNase J. 

Although RNase J deficiency is known to impair 
chloroplast function by allowing the deleterious accumu- 
lation of double-stranded RNA (6), defects in processing 
of functional RNAs could also arise, given that RNase J 
has such a role in bacteria (5,24). Indeed, RNase J-defi- 
cient tissues in both Arabidopsis and N. benthamiana, ac- 
cumulate abnormally sized sense mRNA species as judged 
by RNA gel blot, consistent with a role of RNase J in 
chloroplast RNA maturation (6). Here, we used RACE 
to precisely map the mRNA 5' ends in RNase J-depleted 
tissues of N. benthamiana. 



We initially chose atpH for analysis because its 5' end 
can be correctly generated in vitro by a combination of 
PPR10 and a commercial 5^3' exonuclease (17). As 
shown in Figure 2A and B, 5' end analysis revealed a 
single major band for the wild- type (WT) and the vector 
control (V) tissues corresponding to the previously 
described (13) processed terminus at position —45 (filled 
arrowhead), which abuts the PPR10 binding site. In 
contrast, RT-PCR amplifications of 5' ends from RNase 
J-deficient material revealed additional longer and a few 
shorter termini. Sequence analysis of the longer 5' RACE 
products showed that extended species map between the 
mature WT 5' end of atpH to >400nt upstream, as dis- 
played in Figure 2A. Local secondary structure predic- 
tions suggest that nearly all of the sequenced 5' 
extensions terminate in regions of single-stranded charac- 
ter (Supplementary Figure SI). These findings indicate 
RNase J plays a critical role in the 5' maturation of 
atpH mRNA, perhaps by resecting unstructured 
intergenic regions following endonucleolytic cleavage by 
itself or another enzyme. 

qRT-PCR and RNA gel blots were used for additional 
assessments of sense strand transcript accumulation. The 
qRT-PCR data (Figure 2C) revealed a 2-fold increase of 
sense strand ^//-containing RNAs in RNase J- VIGS 
material compared with control tissues. RNA gel blot 
analysis of the RNase J-deficient tissues (Figure 2D), 
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Figure 1. VIGS-mediated deficiency of RNase J and PPR proteins in N. benthamiana. All plants are shown 18 days post-infiltration. (A) WT, wild- 
type; V, plants infiltrated with the empty pYL170 vector as a negative control; J, RNase J-silenced plants. (B) Plants targeted for PPR proteins 
PPR10, HCF152 and MRL1. (C) VIGS of both RNase J and the indicated PPR protein. 



9144 Nucleic Acids Research, 2013, Vol. 41, No. 19 



100 nt 



-45 



PPR10 

5' extensions: -81,-151,-171,-171,-213, -221,-247, -250, -263, -267, -342, -344, -346, -418, -422 



B 




PPR10 J/PPR10 



Figure 2. Analysis of atpH by 5' RACE, qRT-PCR and RNA gel blot. (A) 15 sequenced atpH 5' termini from RNase J-VIGS plants. The numbers 
shown correspond to nucleotide distances upstream of the atpH translation initiation codon. PPR10, conserved minimal PPR10 binding site (17). 
(B) Gel separation of atpH 5' RACE products. The filled arrowhead indicates the WT atpH 5' end, whereas the open arrowheads designate atpH 
5' extensions that accumulate in RNase J-deficient tissue, as presented in (A). Left lane, 100 bp DNA ladder. WT, wild-type; V, vector control; 
J, RNase J-VIGS. (C) Relative amounts of ^//-containing fragments measured by qRT-PCR. (D) Strand-specific RNA gel blot analysis of atpH. 
An ethidium bromide-stained image of 28 S rRNA is shown below to reflect loading. 



using a strand- specific probe, revealed a background 
smear along with a decrease in the abundance of discrete 
species. One interpretation is that qRT-PCR measure- 
ments reflect a modest increase of atpH mRNAs with het- 
erogeneous 5' ends. The results also indicate that RNase J 
normally limits abundance of atpH mRNA, in the short 
segment of the coding region defined by the qRT-PCR 
primer set. 

To test whether accurate RNase J-mediated atpH 
mRNA maturation requires the cognate PPR protein, 
PPR10, VIGS was used to create PPR10 deficiency with 
or without concomitant RNase J silencing. Targeting of 
PPR 10 (Figure IB) resulted in a chlorotic phenotype, rem- 
iniscent of the maize pprlO null mutant phenotype (13). To 
verify that PPR 10 in tobacco and maize has orthologous 
functions, RNA analysis was performed. Figure 2C shows 
that PPR 10 deficiency does not cause a general loss of 
^//-containing transcripts, whereas Figure 2D clearly 
shows that the levels of transcripts with atpH at their 5' 
ends is strongly reduced, in particular the atpH 
monocistronic mRNA and the atpH-atpF dicistronic 
species (atpF is immediately downstream of atpH). 
Furthermore, because this reduction is not accompanied 
by an increase higher molecular weight bands, we 
conclude that the forms with reduced abundance are 
unstable. PPR 10 deficiency was also associated with a 



shift in migration of the major polycistronic species, 
which may suggest that PPR 10 influences RNA processing 
at secondary sites. The 5' RACE (Figure 2B) showed that 
although the WT 5' end could still be detected, as expected 
from the partial deficiency that VIGS creates, some longer 
5' ends were also present, suggesting that the role of 
PPR10 in defining 5' ends had been compromised. 
Overall, our results agree with findings from maize that 
PPR 10 is necessary for the accumulation of processed 
atpH transcripts (13). 

Simultaneous VIGS silencing was used to investigate the 
relationship between RNase J and PPR 10 during the mat- 
uration of atpH 5' ends, leading to the same whole-plant 
phenotype that was observed following silencing of only 
RNase J (Figure 1). Figure 2B shows that based on 5' 
RACE, dual RNase J/PPR10 deficiency resulted in 
similar extended atpH 5' ends as observed in RNase J- 
depleted plants. However, an effect of PPR 10 deficiency 
was also evident, as ^//-containing transcripts were 
reduced to control levels as judged by qRT-PCR analysis 
(Figure 2C). Apart from transcript abundance, the overall 
similarity of molecular pheno types between RNase J/ 
PPR 10- and RNase J-deficient plants suggests that the re- 
duction of RNase J masks the compromised protective role 
of PPR10. Therefore, RNase J is involved in the 5' matur- 
ation of atpH transcripts bound by PPR 10 in vivo. 



Nucleic Acids Research, 2013, Vol. 41, No. 19 9145 



100 nt 



-393 



psbH 




HCF152 



B 



-393 





psbH-petB/l B -petD - 
petB/l B -petD- 

psbti-petB/D- 

petB/D- 

psbH-petB - 

petB- 




HCF152 J/HCF152 



Figure 3. Analysis of petB mRNA. (A) Schematic of the major bands detected by 5 ; RACE. The filled gray and black arrowheads correspond to the 
similarly-labeled bands in (B). Open arrowheads represent estimated 5 ; ends from products of intermediate sizes. HCF152, conserved minimal 
HCF152 binding site (15). (B) 5' RACE products from RT priming within the petB exon. The filled black arrowhead marks the mature WT 
petB 5' end (25) and the filled gray arrowhead marks products from a precursor transcript that terminates upstream of psbH, corresponding to the 
psbH-petB dicistron. Numbers at right are nucleotide positions relative to the petB translation initiation codon. (C) Relative amounts of 
/?e/i?-containing RNA fragments measured by qRT-PCR. (D) Strand-specific RNA gel blot for the petB exon. I B indicates unspliced petB species. 
An ethidium bromide-stained image of 28S rRNA is shown below to reflect loading. 



RNase J and HCF152 are required for the maturation of 
petB S termini 

To see whether the paradigm verified for RNase J and 
PPR10 holds true in the context of another gene cluster, 
we analyzed the 5' end of petB, a transcript derived from 
the psbB operon, which is protected by the PPR protein 
HCF152 (15). We used the same VIGS strategy as for 
PPR 10, and the expected plant pheno types were 
obtained as shown in Figure 1 — Arabidopsis HCF152 
null mutants are chlorotic and require sucrose for 
growth (14). Figure 3 A and B show that heterogeneous 
and extended petB 5' termini, compared with the WT -44 
position, accumulate in tissue with reduced RNase J. This 
contrasts with the HCF152-depleted tissue, where the 
major species not found in WT (gray arrowhead) corres- 
ponds to the psbH 5' end (—393), indicative of transcripts 
that have not undergone intercistronic processing. The 
presence of residual WT 5' ends is consistent with the in- 
complete silencing typical of VIGS and is in agreement 
with results from the weak mutant allele hcfl52-2 of 
Arabidopsis (14). RNA from RNase J/HCF152-depleted 
material exhibited a 5' end pattern most similar to that 
seen with RNase J deficiency alone, suggesting that as 
for PPR 10, RNase J deficiency can at least partially 
mask HCF152 deficiency. 



Additional assessment of petB mRNAs was performed 
using qRT-PCR and strand-specific RNA gel blot 
analysis. As shown in Figure 3C, RNase J deficiency 
was linked to a 5-fold reduction of ;?<^i?-containing 
RNAs, accompanied by a diffuse RNA gel blot signal 
(Figure 3D). Reduced RNA accumulation is likely to be 
at least in part a secondary consequence of diminished 
PEP activity that is associated with the chloroplast trans- 
lational defects in RNase J-depleted chloroplasts: petB 
transcription is mainly dependent on the upstream psbB 
PEP promoter, in contrast to atpH whose transcription is 
increased in the absence of PEP (26). HCF152-deficient 
plants also exhibited reduced accumulation of ^^-con- 
taining transcripts, including the ones that had been 
cleaved within the psbH-petB intergenic region and espe- 
cially those containing spliced petB. The persistence of 
higher molecular weight petB RNAs in HCF 152- VIGS 
tissue suggests a comparative instability of the spliced 
forms, a phenomenon first noted in Arabidopsis hcfl52 
T-DNA mutants (14). Joint reduction of RNase J and 
HCF 152 led to similar petB mRNA defects as were 
observed for RNase J depletion, where an accumulation 
of 5'-extended petB species (Figure 3B) was accompanied 
by a general decrease in petB RNA abundance (Figure 3C 
and D). These results suggest RNase J and HCF 152 
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operate on the same petB transcripts and are both 
required for proper 5' maturation and stability. 

To see whether there was a differential effect of RNase J 
and/or HCF152 deficiency on unspliced petB species, we 
used 5' RACE initiated by reverse priming within the petB 
intron and RNA gel blots (Supplementary Figure S2). 
Although the major unspliced petB RNAs accumulated 
in HCF152-VIGS material at similar levels as in WT 
tissue, as judged by gel blot, few petB 5' ends at position 
—44 were detected by RACE, with the products being 
dominated by 5' ends upstream of psbH and species 
mapping slightly downstream of the —44 site. In the 
cases of RNase J deficiency, 5' extensions were also 
mapped between —44 and psbH. The effects of HCF152 
deficiency thus differed significantly for spliced versus 
unspliced RNAs, where the —44 5' end was readily 
found for spliced species but barely detectable for 
unspliced RNAs. One interpretation is that spliced 
RNAs have a higher affinity for residual HCF152. The 
small size of the first exon would be consistent with an 
effect of the intron on RNA structure in the HCF152 
binding region (see Figure 3A). Although HCF152 
binding to different regions surrounding, the petB 5' end 
has been assessed (15,27), to our knowledge, a direct com- 
parison of spliced and unspliced substrates has not been 
conducted. 

RNase J and MRL1 are necessary for the processing of 
vbcL 5' ends 

We chose the monocistronic rbcL mRNA as a third model 
to investigate RNase J involvement in 5' end processing, 
as a complement to the polycistronic transcript derivatives 
atpH and petB. Two forms of rbcL mRNA accumulate in 
land plants: a primary transcript and a shorter form 
generated by ribonucleolytic processing (28). These can 
be differentiated before 5' RACE by treatment of RNA 
with TAP because primary transcripts can only be 
amplified following TAP treatment. As shown in 
Figure 4A and B, WT and vector control plants possess 
both primary (—182) and processed (—59) forms of rbcL. 
In the +TAP lanes, the product representing the primary 
form predominates owing to its higher abundance, 
whereas in —TAP lanes, only the processed form is repre- 
sented. Depletion of RNase J led to a strong reduction of 
primary rbcL transcripts, suggesting that PEP-mediated 
transcription of rbcL had diminished, and to the appear- 
ance of 5' extended species relative to the —59 form. The 
species mapping between positions —59 and —182 repre- 
sent incomplete maturation owing to depletion of RNase 
J. The —59 form, but not the primary transcript, is found 
in tobacco where PEP has been genetically eliminated, 
reflecting the presence of additional promoters upstream 
(30), which also accounts for the accumulation of novel 
longer transcripts on the gel blot for RNase J-depleted 
material. 

We next partially silenced expression of the PPR protein 
MRL1, which specifically protects the processed rbcL 
transcript. VIGS material was indistinguishable from the 
WT in terms of plant growth phenotype (Figure IB), con- 
sistent with results from Arabidopsis (16). This is because 



the accumulation of primary rbcL transcripts is unaffected 
by MRL1 deficiency (Figure 4B), while the processed form 
is essentially undetectable. Because of their small size dif- 
ference, the two isoforms of rbcL mRNA did not always 
resolve in RNA gel blot analysis. When RNase J was 
silenced along with MRL1, plants accumulated a ladder 
of rbcL 5' ends that differed little between +TAP and 
—TAP samples, as was the case when only RNase J was 
silenced. The additive effect of reducing both MRL1 and 
RNase J expression was a near-complete loss of —59 rbcL 
termini, as well as disappearance of a discrete species seen 
by gel-blot analysis, that migrates slightly faster than the 
WT transcript (Figure 4C). It is possible that this species is 
defined by a secondary MRL1 binding site. 

To cross-check our conclusions for rbcL, we examined 
psbB 5 f ends. Like rbcL, psbB mRNA features closely 
spaced primary and processed 5' ends, and its primary 
transcript is also PEP-dependent. As shown in 
Supplementary Figure S3, RNase J depletion leads to 
near-absence of the primary transcript and a ladder of 
5' ends ranging from shorter than the processed form to 
well upstream. RNase J, therefore, appears to be required 
for generating fully processed 5' ends from both rbcL and 
psbB primary transcripts, an event originally assumed to 
be catalyzed by a single endonucleolytic cleavage (25). The 
possible trans-factor that protects the processed psbB 
5 r end has not been identified: mutants deficient for the 
Arabidopsis ortholog of Mbbl, which protects psbB 
in Chlamydomonas (31), have RNA defects only for 
psbH (32). 

RNase J processes psbA mRNA derived from 
read-through transcription of an upstream tRNA gene 

Another monocistronic transcript investigated by 5' 
RACE was psbA. Unlike rbcL, psbA has a single predom- 
inant form, assumed to be entirely derived from transcrip- 
tion initiation at its well-studied PEP promoter at —85 
[(33); Figure 5A]. As shown in Figure 5B, RACE 
analysis of TAP-treated RNA from WT and vector 
control plants revealed only the expected primary tran- 
script. When processed species were specifically amplified 
(minus-TAP lanes), smaller products were observed, the 
major ones being 50-75 bp shorter than the primary tran- 
script. Indeed, minor 5' ends at approximately those loca- 
tions were previously noted during mapping of 
Arabidopsis psbA transcripts (35), and shorter fragments 
were also detected during RNA decay assays in isolated 
spinach chloroplasts (36). 

When RNase J was depleted, none of the putative WT 
processed forms were detected, and only a trace amount of 
the primary transcript was found, which appeared to be 
replaced by a product 130 bp longer (Figure 5B). Loss of 
the primary transcript is consistent with reduced PEP 
activity. Similarly, RNA gel blot analysis showed rela- 
tively little of the WT-sized band, along with a 
somewhat less abundant and slightly larger species 
(Figure 5C). The lengthiest of the 5'-extended species 
(filled gray arrowheads in Figure 5A and B) was 
sequenced and found to match precisely the 3' end of 
mature tRNA Lys , which is encoded immediately 
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Figure 4. rbcL 5' end analysis. (A) Graphic representation of rbcL 5' ends detected by RACE. The filled gray and black arrowheads correspond to 
the labeled bands in (B). Open arrowheads represent transcript ends of estimated sizes. MRL1, predicted MRL1 binding site (15). (B) Gel separation 
of 5' RACE products. +, TAP treatment prior to RACE; — , samples not treated with TAP. The filled black arrowhead marks the processed rbcL 
5' end (29) and the filled gray arrowhead designates the primary transcript. Values at right designate nucleotide positions relative to the rbcL 
translation initiation codon. (C) Strand-specific RNA gel blot probed for rbcL. An ethidium bromide-stained image of 28S rRNA is shown below to 
reflect loading. 
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Figure 5. psbA 5' end analysis. (A) Schematic of major psbA 5' RACE 
products. The gray rectangle represents the 3' end of trnK exon 2. The 
filled gray and black arrowheads correspond to the 5 ; ends marked in 
(B). Open arrowheads represent 5' ends of estimated sizes. (B) 5' RACE 
products for psbA. Primary transcript initiating at the PEP promoter 
(filled black arrowhead) (34); read-through transcript from the 
upstream trnK gene (filled gray arrowhead). Positions shown are 
relative to the psbA translation initiation codon. (C) Strand-specific 
RNA gel blot analysis of psbA. The filled gray arrowhead designates 
psbA transcripts arising as byproducts of trnK 3' end maturation. An 
ethidium bromide-stained image of 28S rRNA is shown below to rep- 
resent loading. 



upstream and whose 3' maturation is catalyzed by the 
endonuclease RNase Z (37,38). The fact that amplification 
of these extended transcripts is independent of TAP treat- 
ment confirms their identity as processing products rather 
than derivatives of a second psbA promoter. To estimate 
the frequency of read-through transcription, a run-on 
transcription assay was performed, where lysed chloro- 
plasts were labeled for 5min in the presence of oc- 32 P- 
UTP to extend previously initiated transcripts (39). This 
analysis showed that transcription activity within the psbA 
coding region was ~12 times higher than the activity 
upstream of the PEP promoter in WT chloroplasts 
(Supplementary Figure S4). Assuming similar stabilities 
of transcripts generated through either means, one can 
estimate that 9-10% of psbA transcripts result from 



read-through transcription, and there is no reason to 
believe that those transcripts are not functional. 



DISCUSSION 

The results presented here provide evidence that chloro- 
plast RNase J plays a major role in the 5' maturation of 
mRNA. This could occur by 5'— >-3' exonuclease activity 
following RNase J-catalyzed endonucleolytic cleavage, or 
RNase J could act mainly in a 5'— >-3' mode following 
endonuclease cleavage by other enzymes such as RNase 
E or RNase Z. Figure 6, which builds on a previous model 
(13), illustrates the various modes of RNase J participa- 
tion in chloroplast mRNA 5' end maturation that are con- 
sistent with our data, depending on whether the substrates 
are derived from intercistronic cleavage (A), nearby tran- 
scription initiation (B) or 3' processing of an upstream 
tRNA (C). Figure 6 posits that RNase J acts in partner- 
ship with other endoribonucleases, either RNase E or 
RNase Z. Evidence for RNase Z participation comes 
from results for psbA (Figure 5), whereas results for 
atpH, where RNase J deficiency leads to accumulation 
of RNAs with a range of 5' ends mainly in predicted 
single-stranded regions, implicate the activity of a 
second endoribonuclease (Supplementary Figure SI). 
Concomitant silencing of RNases J and E could in prin- 
ciple test this hypothesis. 

Although the role of RNase J in 5' end maturation is 
readily observed following VIGS, the associated reduction 
of PEP activity must be taken into account. PEP cannot 
be produced without chloroplast translation, and RNase J 
deficiency globally increases antisense RNA accumula- 
tion, which in turn impedes translation (6). Some RNAs 
we examined are more dependent on PEP than others, 
leading to varying effects on their accumulation as 
measured by qRT-PCR or RNA gel blot. Unlike PEP, 
NEP does not depend on chloroplast translation, as it 
lacks chloroplast-encoded subunits, and its activity 
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Figure 6. Model for the involvement of RNase J in the processing of chloroplast mRNA 5' ends defined by PPR proteins. Maturation of a generic 
mRNA (ORF) 5' end is shown here. Precursor transcripts originate from polycistronic (A) and monocistronic (B) transcriptional units, as well as 
read-through transcripts from upstream genes such as tRNAs (C). Processing is initiated by endonucleolytic cleavages by RNase J or RNase E within 
unstructured intergenic regions, or in the case of tRNAs by RNase Z. The resultant 5 ; ends are subsequently trimmed to their mature forms by the 
exonuclease activity of RNase J to PPR protein-bound sites. 



generally increases in the absence of competition from 
PEP (40). It is thus unsurprising that one transcript 
more prevalent in RNase J-deficient material was atpH, 
which is more strongly transcribed when PEP is absent 
(26). As another example, atpB is transcribed by both 
NEP and PEP in WT tissue (41), and in plants, depleted 
of RNase J, the NEP-transcribed atpB persisted while 
PEP-initiated forms were not detected (Supplementary 
Figure S5). 

Mechanism of 5' maturation by RNase J: exonucleolytic 
processing of endonucleolytic products 

A central premise of the model presented in Figure 6 is 
that RNase J behaves as a 5'-*3' exonuclease while 
trimming the 5' ends of mRNA precisely to the borders 
of bound PPR proteins. A recent study revealed that re- 
combinant PPR 10 and a commercial 5 f ^3 f exonuclease 
were sufficient to generate atpH 5' ends that match the 
processed PPRlO-dependent termini found in vivo (17). 
Fine mapping has further shown that the majority of pro- 
cessed chloroplast RNAs in land plants may be generated 
by an analogous mechanism (15). In Chlamydomonas 
chloroplasts, RNA-binding proteins also shield mRNA 
from 5'-*3' exonucleolytic degradation (31,42), suggesting 
the mechanism of 5' stability elements functioning as road- 
blocks against 5'— >-3' exonucleases is well conserved. 
The present work in N. benthamiana shows that RNase J 



is required for the maturation of several mRNA 5' ends 
in vivo (Figures 2-4), which likewise map precisely to the 5' 
borders of predicted PPR protein-binding sites 
(Supplementary Figure S6). 

The exonuclease activity of B.subtilis RNase Jl is 

currently regarded as the enzyme's predominant mode of 
action in vivo (43). B. subtilis possesses an RNA processing 
system related to that of the plant PPR protein-mediated 
maturation, as 5' stability elements can protect kilobases 
of downstream RNA (5). However, Arabidopsis RNase J 
has a relatively robust endonuclease activity in vitro (6). 
This may suggest that in higher plants, the two primary 
endonucleases that initiate mRNA processing are RNase J 
and RNase E, with RNase Z contributing where tRNAs 
are present immediately upstream of protein-coding genes. 
Arabidopsis RNase E null mutants accumulate a few 
lengthy precursor mRNAs, but processed forms are still 
detected, and the abundance of many mRNAs is unper- 
turbed (44,45). These findings support the notion that 
RNase E is not the sole endoribonuclease that processes 
chloroplast mRNAs. Although it is tempting to speculate 
that the 5'-extended mRNA species that accumulate on 
depletion of RNase J arise directly from RNase E cleav- 
ages, we cannot rule out that some proportion of them 
arise due to stalled exonucleolytic processing by residual 
RNase J in VIGS tissue or from activity of yet another 
enzyme. 
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Arguing against the hypothesis that RNase J carries out 
rate-limiting endonuclease cleavages that initiate mRNA 
processing is the fact that lengthy mRNA precursors 
were not found to accumulate in RNase J-VIGS tissue 
(Figures 2D, 3D, 4C and 5C; Supplementary Figures 
S2C and S3C). As VIGS can be assumed to only partially 
deplete plants of RNase J, it is plausible that residual 
RNase J was sufficient to mediate near-normal 
endonucleolytic cleavage. In this respect, recent structural 
analyses of B. sub t His RNase Jl (46) and Thermus 
thermophilus RNase J (47) suggest that complete or 
partial subunit dissociation is necessary for the positioning 
of an internal RNA segment across the enzyme's catalytic 
cleft for endonucleolytic cleavage. If this is also the case 
for the chloroplast enzyme, diluting the concentration of 
RNase J by VIGS would promote greater endonuclease 
activity, potentially allowing proper initiation of mRNA 
processing even while its exonucleolytic activity was 
impaired. 

A prevalent mechanism of chloroplast mRNA 5' end 
maturation involves RNase J and PPR proteins 

RNase J is involved in the maturation of mRNAs with 
diverse origins: atpH and petB are derived from larger 
transcriptional units, rbcL and psbB are primary tran- 
scripts with processed isoforms, and a previously unrec- 
ognized form of psbA results from 3' maturation of an 
upstream tRNA. In addition to its role in processing 
sense transcripts, RNase J degrades transcriptional read- 
through products that are prolific in chloroplasts (6), sug- 
gesting that the enzyme operates on a global scale and 
with little sequence specificity. Gene-specific RNA- 
binding factors, such as PPR proteins, are therefore 
critical in stabilizing functional mRNAs. By precisely 
designating the ends of chloroplast RNAs that are 
destined for maturation, PPR proteins serve as molecular 
barriers against the non-specific degradation by RNase J, 
permitting compatibility of the enzyme's roles in both 
asRNA cleanup and mRNA processing (15). Stem-loops 
may also serve as barriers to RNase J and define 5' ends, 
although we did not demonstrate such a mechanism in this 
study. 

Beyond the PPR proteins investigated here (PPR10, 
HCF152 and MRL1), many other PPR, PPR-like and 
sequence-specific proteins bind mRNA termini and 
confer stability. In Chlamydomonas, MCA1 and MCD1 
protect pet A and petD mRNAs from 5'— ^ exonucleolytic 
decay (42,48). Another chloroplast RNA-binding protein, 
HCF107, serves to define psbH mRNA 5' termini by re- 
sisting 5'— >-3' exonuclease degradation in Arabidopsis and 
Chlamydomonas (31,49,50). Recent deep sequencing of 
small RNAs (sRNAs) has greatly augmented the ability 
to identify transcripts that undergo protein-mediated pro- 
cessing of the type studied here because PPR or analogous 
proteins leave in vivo sRNA footprints (15,51). Not only 
do sRNAs validate the known protective functions of 
PPR 10 and the other PPR proteins studied here but 
conserved sRNAs also map to the 5' ends of psbB and 
psbA (Supplementary Figure S6). This suggests the 



involvement of as-yet uncharacterized RNA-binding 
proteins in 5' end maturation of these transcripts. 

The psbA case is intriguing because this mRNA has 
been presumed to arise uniquely from transcription initi- 
ation and thus harbor a protective triphosphate at its 5' 
end; an assumption supported by the TAP-dependence of 
RT-PCR amplification shown in Figure 5B, and earlier 
data showing that the 5' end could be labeled by a 
capping enzyme (52). Our data, however, suggest that 
some portion of psbA mRNA results from the processing 
of trnK read-through transcripts, and this RNA would 
lack a 5' triphosphate. As all the —85 5' ends are TAP 
dependent, RNase J may convert the downstream 
byproduct of RNase Z-catalyzed trnK 3' maturation to 
the slightly shorter RNAs seen in minus-TAP lanes 
(Figure 5B). If one accepts, however, that the sRNA foot- 
print results from protection of the abundant psbA 
mRNA 5' end by an RNA-binding protein, then the tran- 
scripts produced by RNase J should likewise be protected. 
Isolation and analysis of this predicted psbA binding 
factor may help to resolve this ambiguity. 
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